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Probing the near infrared stellar population of Seyfert galaxies 
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ABSTRACT 

We employ IRTF SpeX NIR (0.8/im-2.4//m) spectra to investigate the stellar population 
(SP), active galactic nuclei (AGN) featureless continuum (FC) and hot dust properties in 9 
Sy 1 and 15 Sy 2 galaxies. Both the STARLIGHT code and the hot dust as an additional base 
element were used for the first time in this spectral range. We found evidence of correla- 
tion among the equivalent widths (W\) Si 1 1.59/xm x Mgl 1.58^m, equally for both kinds 
of activity. Part of the WNai2.2i/«n an d IVco 2.3pm strengths may be related to galaxy incli- 
nation. Our synthesis shows significant differences between Sy 1 and Sy2 galaxies: the hot 
dust component is required to fit the A' -band spectra of ~90% of the Sy 1 galaxies, and only 
of ^25% of the Sy 2; about 50% of the Sy 2 galaxies require a FC component contribution 
>20%, while this fraction increases to about 60% in the Sy 1; also, in about 50% of the Sy2, 
the combined FC and young components contribute with more than 20%, while this occurs 
in 90% of the Syl, suggesting recent star formation in the central region. The central few 
hundred parsecs of our galaxy sample contain a substantial fraction of intermediate-age SPs 
with a mean metallicity near solar. Our SP synthesis confirms that the 1.1 /im CN band can 
be used as a tracer of intermediate-age SPs. The simultaneous fitting of SP, FC and hot dust 
components increased in ~ 150% the number of AGNs with hot dust detected and the mass 
estimated. The NIR emerges as an excellent window to study the stellar population of Sy 1 
galaxies, as opposed to the usually heavily attenuated optical range. Our approach opens a 
new way to investigate and quantify the individual contribution of the three most important 
NIR continuum components observed in AGNs. 
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1 INTRODUCTION 

A key issue in modern astrophysics is to understand the origin of 
the energy source that powers the continuum and line-emitting 
gas in active galactic nuclei (AGN). The current paradigm pro- 
poses that accretion of material onto a supermassive black-hole 
located at the centre of the galaxy is the mechanism responsible 
for the observables associated with the AGN phenomena. In 
addition to this central engine, observational evidence over the 
past decade has shown that massive star forming regions are 



starburst clusters coexist in the nuclear region of galaxies. 
Currently, there is ample evidence indicating that both the 
active nucleus and starburst might be related to gas inflow, 
probably triggered by an axis-asy mmetry perturbation like bars, 
mergers or tidal interactions dShlosman. Frank & Beg elman 



IShlosman. Begelman & Frank 
iKnapen. Shlosman & Peletied 



19901; 
l200Ct 



Maiolino et al 



iFathi et al 



Riff el et al .1 2008a), provid ing support to t h e so- 



an AGN (e.g. |Tokunaga et alj 


1991; iMizutani, Suto & Maiharal 


19941; Imanishi & Dudlev 2000 


; Storchi-Bergmann et al. 2000; 


Imanishi 2002 ; 


Rodriauez-Ardila & Viegasl 120031: Riffel et al.1 


20071; iDorsetal. 


2008). In 


this scenario, black-holes and 



called AGN-starbur s t connection (Norman & Sc ovilld 
Terlevich et al. |ll990l: iHeckman et all ll 9971; ICid Fernandesl 



Gonzalez-Delgado et a l. 1998; Veilleux 2o6b1; iFerrarese & Merrittl 
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2000l : lHeckmarJl2004lRiffel et alj|2008al and references therein). 

Another line of thought, however, claims that this AGN- 
starburst connection could be incidental, as m any Seyferts do 
not sh ow any evidence of starburst activity (e.g. iFilippenko et al .1 
1 1993b , and optical spectroscopic studies of large samples do not 
indicate that star bursts are m ore common in Seyferts than in 
norma l galaxies dPoggelll989h . In addition, ICid Fernande s et al .1 
d2004 hereafter CF04) studied the optical st ellar po p ulatio n (SP) 
of 65 Sy 2 and 14 other galaxies from the Ijoguetl d200ll) sam- 
ple. They concluded that the star formation history (SFH) of 



2 Riffel et al. 



the Sy 2 galaxies is remarkably heterogeneous. These results are 
similar to those obtained in the study of the UV and optical 
SP o f Seyfert g a laxies (mostly Sy 2 ) available in the literature 



presented and discussed in Sect. [5] The final remarks are given in 
Sect.[U 



(e.g. iBicd 119881; IS chmitt 



lostly Sy I ) available in the literature 
et al j 1 19961; ICid Fernandes et all 1 19981; 



Gonzalez-Delgado et all 1998 ; Bonatto et all 200d; Raimann et all 



20031 ; iGonzalez-Delgado et alj|2004 ICid Fernandes et alj|2005bl 2 THE DATA SET 



CF04, and references therein). 

To determine if circumnuclear SPs and nuclear activity are 
closely related phenomena, or if they are only incidental, it is 
of utmost importance the correct characterisation of the former, 
since a substantial fraction of the energy emitted by a galaxy 
in the optical to NIR domain is starlight. Moreover, the anal- 
ysis of the stellar content provides information on critical pro- 
cesses such as the star formation episodes and the evolution- 
ary history of the galaxy. In this respect, the use of NIR fea- 
tures in the study of SP i s not recent, dating back to nearly three 
decades ago. For example. lRieke"et "aH dl980n employed NIR spec- 
troscopy to study NGC 253 and M 82. They report the detection 
of a strong 2.2/xm CO band, suggesting the presence of a dom- 
inant population of red giants and supergiants in the nuclear re- 
gion of both sources. Since then, other authors have used the 
NIR to study star formation, in most cases ba s ed on the 2.2^tm 
CO band (e.g.lOriglia. Moorwood & 01ivalll993l : loiiva et al J 19951 



lEngelbracht et alj 19981 ; Lancon et al.l200l|. and references therein) 



or ph otometric methods (e.g. lMoorwood & G lass 1982; lHunt et al .1 
120031) , 

One reason to use the NIR to study the SP of AGNs is that it 
is the most convenient spectral region accessible to ground-based 
telescopes to probe highly obscured source s. However, tracking th e 
star formation in the NIR is complicated dOriglia & Olivall2000h . 
Except for a few studies such as th ose based on the Br7 emis- 
sion or the CO(2-0) firs t overtone (e.g. lOriglia. Moorwood & Olival 
ll993l ; lQliva etalJl995h . the SP of the inner few hundred parsecs of 
active galaxies in the NIR remains poorly known. Because stel- 
lar absorption features in the NIR are wi dely believed to p rovide 
a means for recognizing red supergiants dOliva et al 1 ll995h . they 
arise as prime indicators for tracing starbursts in galaxies. Besides 
the short-lived red supergiants, the NIR also includes the contri- 
bution of thermally- pulsating asymptotic giant branch (TP-AGB) 
stars, enhanced in y oung to intermediate a ge stellar populations 
(0.2 «: t 2 Gvr. lMarasto"nlll998l . [20051) . The TP-AGB phase 
becomes f ully developed in star s with degenerate carbon oxygen 
cores (see llben & Renzinll ll983;. for a review). Evidence of this 
population in the optical is usually missed, as the most prominent 
spectral features associated with this population fall in the NIR 
dMarastonll2005h . 

With the new generations of Evolutionary Population Synthe- 
sis (EP S) models , which include a proper treatment of the TP-AGB 
phase dMarastor] l2005h . it is now possible to study in more de- 
tail NIR SP of galaxies. According to these models, t he effects 
of TP -AGB stars in the NIR spectra are unavoidable. IMarastonI 



2005h models, by inclu ding empirical spectra of oxygen-rich stars 
Lancon & Woo3l200Cl) . are able to foresee th e presence of NIR 



absorption features such as the 1 . 1 /im CN band jRiffel et alj2007h . 
whose detection can be taken as an unambiguous evidence of a 
young to intermediate age SP. 

Given the above, we feel motivated to carry out the first de- 
tailed study of the SP in active galaxies in the NIR using the entire 
0.8-2.4/im spectral range. This paper is structured as follows: The 
data are presented in Sect. [2] NIR spectral indices are presented 
in Sect. [3] In Sect. [4] we describe the fitting method. Results are 



We chose for this work a subsample of 24 fr om the 47 AGNs with 
NIR spectra published bv lRiffel et all d2006l) . The selected targets 
display prominent absorption lines/bands and are listed in Tab[T] 
All spectra were obtained at the NASA 3 m Infrared Teles cope 
Facility (IRTF). The SpeX spectrograph teavner et al]|2003l) was 
used in the short cross-dispersed mode (SXD, 0.8-2.4/im). The de- 
tector used was a 1024x 1024 ALADDIN 3 InSb array with a spa- 
tial scale of 0.15'7pixel. A 0.8"xl5" slit was employed giving a 
spectral resolution of 360 kms" 1 . The radius of the central inte- 
grated region is few hundred parsecs for almost all sourcefl For 
more details on the instrumental configurati on, data reduction, cal- 
ibration processes and integrated region see IRiffel et all feood) . A 



rapid inspection of Figs. 9, 10 and 12 in lRiffel et alj ( feOOq ) shows 
that all the chosen spectra are dominated by strong abs orption fea- 
tures, th e most prominent ones are identified in Fig . 1 o f iRiffel et alj 
(2008B- 



3 NEAR INFRARED SPECTRAL INDICES : DIRECT 
MEASUREMENTS 

For comparison with published works and future NIR stellar pop- 
ulation studies, we compute the equivalent widths (W\) of the 
NIR absorption features as well as selected continuum fluxes (F\) 
in regions free from emission/absorption, normalized to unity at 
1. 223^tm. W\ and F\ are measured according to the definitions 
of lRiffel et aT] l l2008bl) . The values of W x and F\ are presented in 
Tabs. [T] and [2] respectively. 

The most studied absorption lines in the literature and present 
in our spectra l range are the Calcium triplet line s (Ccffle.g. 
Terlevich et al] Il99d ; iGarcia-Rissmann et alj 120051 ; IVega et all 



20081 and references therein). 



We have four of our objects (Mrk 573, Mrk 10 66, NGC 21 10 
and Mrk 1210) in common with IVega et alj d2008l) . For the first 
3 objects we measure values larger than those reported by them. 
The difference is due to the differen t in dex definitions. W hile 
we use those of |Bica & Alloin I d 19871) and IVega et al] d2008l) use 
ICenarro et al. I d200ll) . The CaT in Mrk 1210 occurs near the de- 
tection lim i t in ou r spectrum, which thus preclude measurements. 
IVega et alj d2008h report a value of 6.22±0.48A for this feature. 
In addition, the diffe rences between our measurements and those 
of lVegaetal.1 d2008l) can be related to the fact that they measure 
WcaT in their synthetic spectra, and therefore our measurements 
represent better the conditions observed in actual galaxy spectra. 
Other explanation for the discrepancies probably lies in the differ- 
ent apertures used in both w orks. They use a slit width of 2.0" 
dGarcia-Rissmann et alj 2005 ) while we use a slit 2.5 times nar- 
rower dRiffel et alj|2006l) . iTerlevich et alj dl99d) report CaT mea- 
surements for NGC 3227, NGC 262 (Mrk 348) and NGC 5953. Our 
WcaT value for the former is consistent with theirs, while for the 
latter two we measure higher values. The discrepancy is again due 



1 Lower than 300 pc for 15 objects, between 300 pc and 500 pc for 6 a nd 
>500 for 3. For more details see col. 10 of Table 1 of lRiffel et all <2Q06|) 

2 Wca,T=WcaTl 0.849/jm + WcaT2 0.854(jm + WcaT3 0.867/*™ 
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Table 1. Equivalent widths measured in the galaxy sample (in A). 



Object/Ion 


CaT, 


CaT 2 


CiiT 3 


CN 


Al 1 


Nal 


Sil 


Mg l 


Sil 


CO 


Nal 


Cal 


CO 


CO 


CO 


BP b (/xm) 


0.8476 


0.8520 


0.8640 


1.0780 


1.1200 


1.1335 


1.2025 


1.5720 


1.5870 


1.6110 


2.1936 


2.2570 


2.2860 


2.3150 


2.3420 


BP r (^m) 


0.8520 


0.8564 


0.8700 


1.1120 


1.1300 


1.1455 


1.2200 


1.5830 


1.5940 


1.6285 


2.2150 


2.2740 


2.3100 


2.3360 


2.3670 


center(/xm) 


0.8498 


0.8542 


0.8670 


1.0950 


1.1250 


1.1395 


1.2112 


1.5771 


1.5894 


1.6175 


2.2063 


2.2655 


2.2980 


2.3255 


2.3545 


to 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


(14) 


(15) 


(16) 



Seyfert 2 



NGC 262(f) 
Mrk 993 
NGC 591 
Mrk 573 
NGC 1144 
Mrk 1066 
NGC 1275 
NGC2110(t) 
ESO428-G014 
Mrk 1210 
NGC 5728 
NGC 5929 
NGC5953({) 
NGC 7674 
NGC 7682 



3.82±0.44 5.92±0.44 4.11±0.42 

1.41±0.15 3.46±0.13 2.73±0.12 - - 2.57 ±0.01 

5.80±0.39 4.57±0.29 4.72±0.17 16.95±0.13 

3.80±0.18 4.04±0.18 3.08±0.17 14.05±0.15 

4.72±0.43 5.88±0.35 4.28±0.21 18.90±0.21 2.14±0.01 3.41 ±0.01 

4.05±0.11 4.50±0.08 3.61±0.07 13.82±0.14 - 1.90 ±0.01 



5.28±0.28 7.47±0.90 4.59±0.20 20.00±0.20 
4.31±0.06 6.30±0.04 3.85±0.01 



2.28 ±0.01 
2.87 ±0.07 



3.56±0.19 5.58±0.17 4.06±0.15 15.42±0.25 - 1.00 ±0.91 

4.68±0.23 6.50±0.17 5.26±0.13 12.65±0.35 1.95±0.01 1.98 ±0.01 

2.24±0.51 3.42±0.43 3.57±0.31 - - 3.20 ±0.36 

3.22±0.20 4.53±0.17 1.49±0.30 11.00±0.23 



1.41±0.27 
3.58±0.13 
3.83±0.07 
2.86±0.01 3.92±0.01 
3.30±0.01 4.27±0.04 
4.58±0.01 

3.78±0.05 
4.77±0.01 
6.54±0.19 
5.41±0.17 
2.17±0.91 4.27±0.49 
1.62±0.01 3.98±0.08 
2.35±0.14 
4.04±0.08 



2.17±0.05 
2.22±0.05 
2.47±0.02 
2.06±0.01 
2.31±0.13 
2.25 ±0.02 

1.81±0.0I 
2.57±0.01 
2.82±0.16 
3.74±0.29 
1.70±0.19 
1.90±0.04 
1.98±0.03 
2.92±0.02 



2.43±0.04 
2.81 ±0.03 
2.76±0.05 
2.81±0.01 
3.48±0.13 
2.71±0.10 
3.20±0.28 
3.40±0.03 
3.35±0.13 
5.73±0.12 
3.52±0.21 
3.61±0.15 

2.70±0.02 



1.04±0.04 
3. 03 ±0.04 
4.51±0.2I 
2.77±0.04 
4.01 ±0.02 
3.68±0.01 
3. 93 ±0.08 
3.30±0.04 
4.80±0.06 

8.74±0.42 
5.43±0.17 
3.82±0.03 



2.75±0.01 0.74±0.03 1.78±0.02 

9.04±0.01 5.24±0.03 8.34±0.03 

5.39±0.13 9.83±0.11 1.39±0.03 10.14±0.09 

1.48±0.03 7.87±0.43 - 4.26±0.13 

1.70±0.09 10.63±0.09 5.76±0.01 

3.11±0.02 9.98±0.28 4.81±0.13 7.87±0.01 

2.49±0.07 

0.75±0.01 4.78±0.05 2.32±0.10 4.20±0.04 

2.31±0.10 12.34±0.34 5.10±0.25 I1.80±0.11 

1.07±0.16 7.73±0.15 - 3.41±0.47 

5.44±0.14 8.I6±0.01 8.98±0.01 10.27±0.01 

3.67±0.36 13.77±0.06 6.78±0.03 I1.42±0.04 

2.67±0.04 13.37±0.20 7.74±0.03 

9.42±0.26 8.33±0.08 



Seyfert 1 



Mrk 334 
NGC 1097 
MCG-5-13-17 
Mrk 124 
NGC 3227 
NGC 4051 
Mrk 291 
Arp 102B 
Mrk 896 



1.12±0.11 3.40±0.06 1.78±0.08 3.13±0.09 1.66±0.04 I.85±0.07 6.14±0.01 2.96±0.00 3.54±0.05 



4.48±0.I7 5.84±0.15 3.29±0.13 
2.35±0.10 3.18±0.O8 3.84±0.08 



1.67±0.06 2.09±0.05 



5.88±0.26 1.79±0.22 



2.40±0.07 4.45±0.()9 



1.46±0.13 1.81±0.13 4.51±0.15 2.62±0.20 3.69±0.10 2.99±0.01 2.04±0.05 

1.91±0.21 - 2.87±0.04 2.25±0.05 3.07±0.06 1.43±0.05 1.97±0.02 

2.19±0.06 1.29±0.30 2.12±0.05 

5.78±0.15 - 3.43±0.03 1.58±0.03 

3.22±0.01 1.93±0.03 3.19±0.08 

2.56±0.01 - 3.43±0.02 - 3.80±0.11 

2.12±0.39 1.38±0.07 3.82±0.05 1.97±0.03 2.14±0.02 



2.36±0.03 1.35 ±0.04 
1.13±0.01 0.73±0.01 



9.45±0.20 
7.95±0.11 

6.24±0.01 
3.69±0.11 
11.90±0.05 
4.61±0.05 



5.16±0.12 10.32±0.02 
3.87±0.03 6.55±0.05 



2.79±0.00 
0.66±0.08 



4.47±0.03 
3.50±0.02 



1.69±0.14 1.24 ±0.25 2.I6±0.01 0.67±0.08 0.53±0.01 1.23±0.06 l.89±0.05 



Table Notes: BP^ and BP r are the blue and red bandpass boundaries. (J) The continuum of Wc a T is affected by spurious emission. If we use only the points free from emission to set the continuum we get 
3.28±0.20, 5.03±0.18 and 3.36±0.17 for CaTi, CaT 2 , CaT 3 , respectively. 



to the different contin uum definitions^ Note that, for NGC 3227, 
iTerlevich etafl dl990t) used an optional blue continuum band at A 
8582A, wich provides a continuum slope very similar in both stud- 
ies, wich thus explains the similarity in WcaT- 

In Fig. Q] we compare VKivigi 1.58^im with 
W^sn 1.59/im. These two absorption lines are corre- 
lated as W / sii=(0.53±0.02)W / Mg i, with CC=0.67 (or 
Wsii=(0.71±0.21)W^°gi ±0 ' 2) , CC=0.69). The correlation 
between Wx of these two lines suggests that almost all the objects 
studied here follow the same chemical enrichment. This hypothesis 
can be associated with the fact that Si and Mg are more abun- 
dant than the other a elements in the Galactic Globular Cluster 
NGC 6 121, which is located near the Galaxy centre (7=350.97°, b= 
15.97°. llvIariano et alj|2008h . 

We have tried diagrams involving other NIR absorption lines 
(or bands), but o nly weak correlatio ns are found. Since the Fx val- 
ues are studied in lRiffel et all d2006h we only present them in Tab. [2] 
for comparison purposes. 

Our sample is composed mostly of spiral galaxies (Tab. 0, 
which tend to in crease the optical WnuI with the inclination b/a 
jBica etalj|l99lh . This could be the case of NGC 5728, which has 
a low b/a, high Vt^Nai and displays the highest y4 v =3.09 mag value 

3 In order to compare our WcuT values with those of ITerlevich et alj 
Jl990l) we compute the WcuT of our NGC 5953 spectrum using their 
continuum definitio. The values thus obtained are: WcaTl = 1.7 A 
W CaT2 =43A and W Ca y3=3.8A. which are very similar to those of 
ITerlevich et"aT1 ll990l) 



of our sample. In Fig.|2]we investigate the relation between galaxy 
inclination and Wx of NIR absorption lines. The correlation sug- 
gests that the interstellar medium plays an important role in the NIR 
VKNafl O n me other hand, Wco is weakly correlated with galaxy 
inclination (Fig. [2] top right), suggesting a weaker dependence on 
the interstellar medium. No correlation of WcaT with b/a is found. 
In addition, the correlation between Nal and Wco 2.3/^m with b/a 
has no evident relation with Seyfert type (Fig. [2}. 



4 SPECTRAL SYNTHESIS 

In this section we study the NIR SP of our galaxy sample, fitting 
the underlying continuum of the 24 AGNs in the spectral range 
between 0.8/itn and 2.4/itn. 

4.1 The base set 

Clearly, the most important ingredient in the SP synthesis is the 
spectral base set, bj.x- An ideal base of elements should cover the 
range of spectral properties observed in the galaxy sample, provid- 
ing enough resolution in ag e and metallicity to p roperly address the 
desired scientific question l lSchmidt et aljl99~lt CF05). 

4 The better NIR correlation can also be associated with the fact that the 
NIR probes SPs deeper in the dust and therefore the contribution of the Na I 
is enhanced by the intrinsic stellar light. 
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Table 2. Continuum fluxes, normalized to unity at 1.223/im. 



Galaxy 





0.81 


0.88 


0.99 


1.06 


1.22 


1.52 


1.70 


2.09 


2.19 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


Seyfert 2 


NGC262 


128 


1.12 


1.11 


1.01 


1.00 


0.94 


0.95 


0.95 


0.92 


Mrk 993 


1.25 


1.26 


US 


1.13 


1.00 


0.80 


0.71 


0.71 


0.37 


NGC591 


1.08 


1.05 


1.05 


1.01 


1.00 


0.76 


0.67 


0.67 


0.35 


Mrk573 


1.26 


1.27 


1.23 


1.15 


1.00 


0.79 


0.72 


0.72 


0.-15 


NGC1144 


1.25 


1.25 


1.21 


1.16 


1.00 


0.76 


0.67 


0.67 


0.34 


Mrk 1066 


1.05 


1.1 1 


1.10 


1.08 


1.00 


0.83 


0.80 


0.80 


0.49 


NGC 1275 


1.47 


1.29 


1.20 


1.16 


1.00 


0.87 


0.81 


0.73 


0.73 


NGC2110 


0.97 


0.95 


0.9S 


1 .05 


1.00 


0.93 


0.89 


Q.89 


0.65 


ESO428-G014 


1.18 


1.17 


1.12 


1.08 


1.00 


0.79 


0.72 


0.72 


0.37 


Mrk 1210 




1.25 


1.27 


1.16 


1.00 


0.85 


0.82 


Q.82 


0.65 


NGC 5728 






1.01 


1.00 


1.00 


0.82 


0.76 


0.76 


0.42 


NGC 5929 


1.25 


1.24 


1.15 


1.10 


1.00 


0.80 


0.73 


0.73 


0.35 


NGC 5953 


1.28 


1.26 


1.17 


1.12 


1.00 


0.76 


0.69 


0.69 


0.33 


NGC 7674 


1.21 


1.12 


1.08 


1.02 


1.00 


0.97 


1.02 


1.02 


1.10 


NGC 7682 


1.04 


I.Od 


1.07 


1.05 


1.00 


0.83 


0.74 


0.74 


0.35 


Seyfert 1 


Mrk 334 


1.07 


1.07 


1.0S 


1 .06 


1.00 


0.93 


0.91 


0.91 


0.69 


NGC 1097 


1.08 


1.25 


1.16 


1.15 


1.00 


0.81 


0.74 


0.74 


0.39 
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1.43 


1.41 


1.25 


1.18 


1.00 


0.79 


0.73 


0.73 


0.44 


Mrk 124 
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1.00 


1.00 


1.07 


1.06 


1.06 


1.00 


NGC 3227 


1.28 


1.25 


1.16 


1.08 


1.00 


0.88 


0.84 


0.84 


0.58 


NGC 4051 




1.27 


1.15 


1.07 


1.00 


0.89 


0.86 


Q.86 


0.73 


Mrk 291 


1.73 


1.51 


1.30 


1.26 


1.00 


0.80 


0.70 


0.70 


0.41 


Arp 102B 


1.48 


1.33 


1.24 


1.23 


1.00 


0.77 


0.69 


0.69 


0.40 


Mrk 896 


1.31 


1.20 


1.12 


1.07 


1.00 


0.95 


0.93 


0.93 


0.73 
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Figure 1. Diagram of WMgl 1.58/im versus Wsn 1.59/^m. The dashed 
line is the a liner correlation. Open symbols are Sy 1 and filled Sy 2. 



Table Notes: The errors on F\ are 3% in all cases. 



One improvement here over previous approaches that at- 
tempted to describe the stellar content of active galaxies using NIR 
spectroscopy is the inclusion of EPS models that take into account 
the eff ects of TP-AGB stars. Accordingly, we use as base set the 
EPS of lMarastorJ J2005l> . The SSPs used in this work cover 12 ages, 
t= 0.01, 0.03, 0.05, 0.1, 0.2, 0.5, 0.7, 1, 2, 5, 9, 13 Gyr, and 4 metal- 
licities, namely: Z= 0.02 Zq, 0.5 Zq, 1 Zq and 2 Zq, summing up 
48 SSPs. 

When trying to describe the continuum observed in AGNs, 
the signature of the central engine cannot be ignored. Usually, this 
component is represented by a featureless continuum (FC, e.g. 
lKoskilll978l CF04) of power-law form that follows the expres- 
sion F v oc v . Therefore, this component was also added to 
the base of elements. The contribution of this continuum (in per- 
centage) to the flux at Ao (1.223/im) is denoted by FC in Tab.[3| 
According to th e unified model (e.g. lAntonucci & Millerlll985l : 
I Antonuccil 1 1 9931) . the FC in Sy 2 galaxies (if present) is due to 
scattered light from the hidden Seyfert 1 (Sy 1) nucleus. How- 
ever, the reader must bear in mind that a common problem in the 
study of the SP of Seyfert galaxies (especially in the optical) is that 
a reddened young starburst (t g;5My r) is indistinguishable from 
an AGN-type continuum (see Sec.|572"llCid Fernandes & Terlevichl 
ll995l : IStorchi-Bergmann et alJl2000L CF04). To avoid this problem 
we have not included very young SSPs in our base (this point is 
also discussed in Sec. |5.2t . 

In the spectral region studied here, hot dust plays an impor- 
tant role in the continuum em ission of active galaxies. Previous 
studies (i,e,. lRiffel et al.ll2006l for instance) report a minimum in 
the continuum emission around 1.2/im, probably associated with 
the red end of the optical continuum related to the central engine 



and the onse t of the emission due to reprocessed nuclear radia- 
tion by dust tearvainisll 19871: iThompsoril 19951: iRudv et alj[200ol: 
Rodriguez- Ardila & Viegaj 120031 : iRodriguez-Ardila & Mazzalavl 
120061 : iRiffel et al.l2006l) . In order to properly account for this com- 
ponent, we have included in our spectral base 8 Planck distributions 
(black-body-BB), with T ranging from 700 to 1400 K, in steps of 
100 K. The lower limit in T is due to the fact that lower tempera- 
tures are hard to detect. Even a small fractional contribution would 
require a sizeable amount of dust. In order to illustrate this point, 
we plot in Fig. [3] combinations of different BB distributions with 
the synthetic template of the starburst galaxy NGC 7714, derived 
by IRiffel et ail d2008bl) .The combination was made summing up, 
in the whole spectral range, increasing fractions of the dust com- 
ponent from 0% to 100%. Thus, we start with the pure NGC 7714 
synthetic spectrum and end with a pure BB distribution according 
to: 

f =( i -i4)^ + (i4) Fs - (i) 

where / is the percentual flux, which we vary in steps of 1%; 
Fsp is the flux of the synthetic spectrum of NGC 77 14 normal- 
ized to unity at 1.223/im and Feb is the BB flux also normalized 
at the same wavelength. As can be observed in Fig. [3] small frac- 
tional contributions of cool (<700K) dust can significantly alter 
the strength of the absorption lines of the K band spectrum, and 
therefore, are very hard to be detected in our spectral range. Hotter 
BB distributions are not used because 1400 K is very close to the 
sublimation tempe rature of graphi t e grains (likely the main con- 
stitue nt of the dust.lBarvainis 1987l: lRodriguez-Ardila & Mazzalavl 
l2006h . 
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4.2 The method 

The second most important ingredient of a stellar population syn- 
thesis is is the code that will suitably combine the individual com- 
ponents of the base elements to construct the final model that will 
represent the observed continuum. Here, as a synthesis code, we 
use for th e first time - in this spec t ral range - the STARLIGHT 
software fcid Fernandes et al.l 120041 ]2005af iMateus et all 120061; 
I Asari et all 2007 ; Cid Fernandes et alfeoOlj ). This code is well de- 
scribed in lCid Fernandes et al.l (CF04. l2005al hereafter CF05). In 
summary, STARLIGHT mixes computational techniques originally 
developed for semi empirical population synthesis with ingredi- 
ents from evolutionary synthesis models(CF05). Basically, the code 
fits an observed spectum Ox with a combination, in different pro- 
portion s, of N+ s i ngle s tellar populations (SSPs). Due to the fact 
that the lMarastonl j2005l) models includ e the effect of the TP-AGB 
phase, crucial to model NIR SP (see iRiffel etal] l200l l2008tl; 
?), we used this EPS models as the base set for STARLIGHTQ 
Extinction is modeled by STARLIGHT as due to foreground dust, 
and parametrised by the V-band e xtinction Ay. We use the CCM 
dCardelli. Clayton & Mathisll 19891) extinction law. Essentially, the 
code solves the following equation for a model spectrum M\ 
(CF05): 



M A = M, 



(2) 



were bj,\ r\ is the reddened spectrum of the jth SSP normalized at 
A ; r\ = 10- 4(Aa_Aao) is the reddening term; M\ is the syn- 
thetic flux at the normalisation wavelength; x is the population vec- 
tor; ® denotes the convolution operator and G(v*, ct*) is the gaus- 
sian distribution used to model the line-of-sight stellar motions, it 
is centred at velocity v* with dispersion cr*. 

The fit is carried out with a simulated annealing plus Metropo- 
lis scheme, which searches for the minimum of the equation 
(CF05): 



A 



(3) 



where emission lines and spurious features are masked out by fixing 
w x =0. For more details on STARLIGHT see CF04 and CF05. 



5 SYNTHESIS RESULTS AND DISCUSSION 

We present in Figs. [4] to [9] the results of the spectral synthesis fit- 
ting procedure. For each galaxy the top panel shows the observed 
and synthetic spectra normalized to unity at 1.223^m. Note that in 
all cases the synthetic spectrum was shifted by a constant for visu- 
alisation purposes. The bottom panel shows the residual spectrum 
Ox — M\. As expected, the residual is dominated by the nebular 
emission. The analysis of the emission lines free from the SP con- 
tamination is beyond the scope of this paper and is left for a forth- 
coming work (Riffel et al., 2009 in preparation). The results of the 
synthesis are summarised in Tab. [3] The quality of the fits are mea- 
sured by the \ 2 (column 13jf|and the adev (column 14) parameters. 
The latter gives the percentage mean deviation \0\ — M\\/0\ over 
all fitted pixels. 



5 As default base set, the code use the SSPs of lBruzual & CharloJ fe003l) . 

6 Which, in fact, is the \ 2 divided by the number of A's used in the 
fit. Reliable fits are obtained when x 2 ~1 (CF04). for more details see 
STARLIGHT manual available at http://www.starlight.ufsc.br 
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Table 3. Synthesis Results. For more details see text. 



Galaxy 


b/a* 


Morphology^ 


FC 


BB C 


BB h 


Xy 


XI 


X.Q 


my 


mi 


mo 
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adev 


Av 


(log U) L 


(log U) M 


{Z*)l 


<2*>M 








(%) 


(%) 


(%) 


m 


(%) 


(%) 


(%) 


(%) 


(%) 






(mag) 


(yr) 


(yr) 


(*) 


(*) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


(14) 


(15) 


(16) 


(17) 


(18) 


(19) 


Seyfert 2 


NGC 262 


1.0 


SO-a 


20.0 


1.8 


3.6 


0.0 


13.0 


61.2 


0.0 


2.9 


97.1 


1.52 


1.06 


1.37 


9.68 


9.89 


0.011 


0.005 


Mrk993 


0.32 


Sa 


4.7 


0.0 


0.0 


0.0 


27.8 


66.8 


0.0 


6.4 


93.6 


0.01 


1.90 


1.32 


9.54 


9.83 


0.013 


0.005 


NGC 591 


0.77 


SO-a 


0.0 


0.0 


0.0 


14.9 


27.4 


57.5 


1.4 


7.5 


91.1 


0.01 


2.87 


1.97 


9.28 


9.84 


0.014 


0.003 


Mrk573 


1.0 


SO-a 


22.4 


0.0 


0.0 


11.9 


52.6 


11.8 


3.4 


45.1 


51.5 


2.98 


1.02 


0.99 


8.93 


9.56 


0.028 


0.024 


NGC 1 144 


0.64 


E 


0.0 


0.0 


0.0 


0.0 


71.3 


26.1 


0.0 


44.9 


55.1 


0.01 


2.22 


1.44 


8.97 


9.32 


0.021 


0.015 


Mrk 1066 


0.59 


SO-a 


17.9 


0.0 


0.0 


4.6 


50.7 


26.8 


0.8 


29.8 


69.5 


2.36 


1.11 


1.54 


9.19 


9.67 


0.025 


0.015 


NGC 1275 


0.77 


SO 


65.8 


0.3 


0.0 


2.8 


0.0 


32.1 


0.1 


0.0 


99.9 


1.15 


0.80 


0.22 


9.68 


9.94 


0.003 


0.001 


NGC 21 10 


0.76 


E-SO 


33.6 


0.0 


0.9 


9.5 


36.3 


18.9 


1.7 


26.0 


72.3 


1.42 


0.75 


1.98 


8.98 


9.70 


0.028 


0.010 


ESO428-G014 


0.58 


SO 


0.0 


0.0 


0.0 


9.2 


49.1 


39.6 


1.1 


15.9 


83.0 


0.03 


2.05 


1.52 


9.13 


9.76 


0.021 


0.005 


Mrk 1210 


1.0 


s? 


12.0 


0.5 


0.0 


9.2 


48.8 


29.1 


0.3 


30.7 


69.0 


2.29 


1.05 


1.14 


9.17 


9.77 


0.013 


0.007 


NGC 5728 


0.58 


Sa 


0.0 


0.0 


0.0 


0.0 


69.5 


27.8 


0.0 


24.4 


75.6 


0.01 


2.89 


3.09 


8.64 


9.51 


0.037 


0.037 


NGC 5929 


0.78 


Sa 


0.0 


0.0 


0.0 


7.3 


24.5 


69.7 


1.4 


11.0 


87.7 


0.01 


2.38 


1.48 


9.37 


9.61 


0.020 


0.022 


NGC 5953 


0.81 


SO-a 


0.0 


0.0 


0.0 


0.0 


33.5 


65.2 


0.0 


7.0 


93.0 


0.01 


2.00 


0.81 


9.46 


9.80 


0.031 


0.023 


NGC 7674 


0.91 


SBbc 


29.2 


2.2 


5.7 


3.4 


24.6 


35.6 


0.7 


15.8 


83.4 


0.03 


1.37 


1.03 


9.28 


9.69 


0.018 


0.010 


NGC 7682 


0.92 


Sab 


0.0 


0.0 


0.0 


0.0 


34.7 


62.2 


0.0 


6.4 


93.6 


0.02 


3.76 


1.89 


9.53 


9.91 


0.020 


0.011 



Seyfert 1 



Mrk 334 


0.70 


Sbc 


23.7 


0.0 


8.2 


12.2 24.1 


31.7 


2.6 


9.9 


87.5 


0.01 


1.45 


1.36 


9.03 


9.71 


0.018 


0.007 


NGC 1097 


0.67 


SBb 


4.3 


0.0 


0.0 


19.8 46.4 


29.8 


2.9 


28.6 


68.5 


1.70 


1.33 


1.29 


8.88 


9.49 


0.024 


0.017 


MCG-5-13-17 


0.67 


E-SO 


6.6 


0.0 


3.9 


22.9 45.7 


21.6 


5.2 


25.1 


69.7 


1.68 


0.81 


0.88 


8.55 


9.47 


0.023 


0.012 


Mrk 124 


0.67 


S? 


36.5 


0.0 


22.6 


3.6 18.3 


19.0 


1.4 


14.8 


83.8 


0.10 


1.29 


0.81 


9.10 


9.73 


0.008 


0.002 


NGC 3227 


0.66 


SABa 


31.0 


0.0 


2.1 


38.4 28.1 


0.0 


10.8 


89.2 


0.0 


0.05 


1.24 


0.99 


8.00 


9.07 


0.011 


0.011 


NGC 4051 


0.75 


SABb 


41.0 


0.0 


6.8 


37.6 14.4 


0.0 


46.7 


53.1 


0.2 


0.03 


1.28 


0.57 


7.53 


8.02 


0.020 


0.010 


Mrk 291 


0.67 


SBa 


6.5 


0.0 


4.5 


0.0 6.1 


85.3 


0.0 


1.1 


98.9 


1.01 


2.06 


0.45 


9.85 


9.94 


0.006 


0.004 


Arpl02B 


0.78 


7 


7.9 


0.0 


3.0 


11.7 35.2 


41.5 


1.0 


12.7 


86.3 


1.61 


1.14 


0.68 


9.08 


9.80 


0.015 


0.005 


Mrk 896 


0.72 


Sab 


27.0 


0.0 


13.7 


4.5 2.3 


53.0 


0.5 


0.6 


98.9 


0.01 


0.90 


0.55 


9.72 


9.96 


0.009 


0.012 



Table notes: ($) From NED. (^) From HyperLeda - Database for physics of galaxies (http://leda.univ- lyonl.fr| IPaturel et al]|2003l) . (*) Abundance by mass 
with Z© = 0.02. 



To take into account noise effects that dump small differ- 
ences between similar spectral components, we followed CF05 and 
present our results using a condensed population vector, which is 
obtained by binning the x into young, xy (tj ^ 5 x 10 7 yr); 
intermediate-age, xi (1 x 10 s ^ tj ^ 2 x 10 9 yr) and old, 
xo (tj > 2 x 10 9 yr) components, using the flux contributions. 
The same bins were used to represent the mass components of 
the population vector; my, mi and mo, respectively. The con- 
densed population vectors are presented in columns 7 to 12 of 
Tab. [3] For more details on vector definition see CF05. We have 
also binned the black-body contributions into two components. The 
cool (BB C ) is obtained by summing up the BB contributions with 
T < 1000 K, and the hot one (BB h ) with T> 1 100 K. These com- 
ponents were defined based on the sublimation te mperatures of sil- 
icate (~1000K) and grap hite (~ 1200 K) grains dBarvainisll 19871 : 
iGranato & Daneselll994h . The condensed black-body vectors for 
each source are listed in columns 5 and 6 of Tab. [3] 

The spectral synthesis shows that the NIR continuum of active 
galaxies can be explained in terms of at least three components: a 
non-thermal continuum, the dust emission and the SP of the cir- 
cumnuclear region. As can be seen in Table [3] and in Fig. [TO] the 
contribution of the latter to the nuclear continuum is higher than 
50% in most objects. Therefore, its study is a critical step in the 
analysis of the continuum emission of Seyfert galaxies. Moreover, 



our results are consistent with the predictions of the unified model 
for AGNs, as the non-thermal continuum and the hot dust emission 
are present in all Sy 1 sources and only in a small fraction of the 
Sy 2s (see also Sec. [33}. In the following sections we provide a 
detailed description of each of these three components in the light 
of the results obtained here and compare them with those obtained 
in other wavelength regions, mostly in the optical. 

5.1 The stellar population 

Binning the population vectors into six component^! left us with 
a coarser but more powerful description of the SFH of our galaxy 
sample. To better quantify the NIR SFH we plot a histogram with 
the flux-weighted and mass-weighted condensed population vec- 
tors in Fig.fjjJ Overall, the NIR SPs are heterogeneous, as in most 
sources the three components contribute significantly to the inte- 
grated flux. However, the contribution of xy is very small (<10%) 
for most of our sample. The intermediate-age component is well 
distributed, with a maximum centred at ~40%. The xo contribu- 
tion is very similar to that of the intermediate age. 

Regarding the mass-weighted components, as expected from 
the results seen above, the contribution of m y is very small 

7 With 3 components representing the star formation episodes. 
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Figure 4. Spectral fits. Each panel shows: (i) at top the flux of the observed spectrum, normalized at unity at 1.223^m, and the synthetic spectra (shifted down 
for clarity); (ii) at bottom the 0\ — residual spectrum. 



(near zero) as can be observed in the right side of Fig. [TT] The 
intermediate-age mass contributions are distributed over all frac- 
tions (from to 100%) but tend to be biased to values lower than 
20%. In contrast, the old component of the mass- weighted vector 
is biased to values higher than 70%. 



According to CF05, if one would characterise the SP mixture 
of a galaxy by a single parameter, it is the mean stellar age. They 
defined it in two ways: the first is weighted by light fraction, 



(logt*) L = ^x j logt j , (4) 

3=1 

and the second, weighted by the stellar mass, 

N, 

<logf*) M = mjlogtj. (5) 

3=1 

Note that both definitions are limited by the age range used in our 
elements base (Sect. [4~T| l and obviously, the FC and BB compo- 
nents were excluded from the sum. The mean stellar ages derived 
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Figure 5. Same as Fig. [4] 



with both definitions are presented in columns 16 and 17 of Tab. [5] 
respectively. 



A secondary parameter to describe the mixed SP is the metal- 
licity. CF05 also defined the light-weighted mean metallicity by 



To better quantify the NIR mean ages we show at the left 
side of Fig. [72l histograms for (logt*)_L and (logt*)j\/. The light- 
weighted mean age of our galaxy sample is biased to an interme- 
diate/old age SP, while for the mass-weighted mean age we clearly 
observe that the old population dominates. As stated by CF05 the 
mass-weighted mean age is a more physical parameter, but it has a 
much less direct relation with the observables. They associated this 
discrepancy with the non-constant stellar M/L ratio. 



j=i 



(6) 



as well as the mass-weighted mean metallicity, which is repre- 
sented by: 



j'=i 



(7) 
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Figure 6. Same as Fig. [4] 



Both definitions are bounded by the -^Zq-IZq range. The light- 
and mass-weighted mean metallicities, estimated for our galaxies, 
are presented in columns 18 and 19 of Tab. [3] respectively. 

We present in the right side of Fig. [12] a histogram for the 
light- and mass-weighted metallicities of our sample. Our results 
point to a mean metallicity solar to above solar, if we consider the 
light-weighted values, while for the mass-weighted mean metal- 
licity our results indicate a sub-solar value. We associate this dis- 
crepancy with the well known age-metallicity degeneracy, i.e. for 
a fixed mass, a high-metallicity SP looks cooler - and older - than 
a low-metallicity SP, thus resulting in a higher M/L ratio. More- 



over, this is consistent with a galaxy chemical enrichment scenario 
in which the young population is enriched by the evolution of the 
early massive stars. In this context, the light-weighted metallicity is 
more sensitive to the young component, while the mass-weighted 
metallicity to the old stellar population. 



5.1.1 Previous studies 

A significant fraction (12 out of 24) of our objects have been sub- 
ject of previous SP studies in the optical and UV regions. In this 
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Figure 7. Same as Fig. [4] 



section we will compare our NIR results with those available in the 
literature. 



• NGC 262: Accordingly to iGarcia- Vargas et alj dl989h most 
of the optical and NIR flux comes from the galaxy bulge SR 
The emission at shorter waveleng ths can be fitted by a power 
law. iGonzalez-Delgad o et alj feOOlh studied the optical SP of this 
galaxy and found a dominant old population. By fitting the nuclear 
continuum with an off-nuclear spectrum of the galaxy bulge as tem- 
plate, they conclude that the contribution of a power law is unneces- 
sary. In the case of an elliptical as a template spectrum, they needed 



to include a contribution ^25% of a power law (Fx <x A" 1 ). They 
argue that such small dilution can be related to a change in the SP of 
NGC 262 with respe ct to the elliptical gala xy template. Similar re- 
sults are obtained bv lRaimann et al.l ( l2003l) . They found 60%, 20%, 
12% and 8% for the old, intermediate, young and FC components, 
respectively As can be observed in Tab. [3] our results agree with 
the optical studies. 

8 Note that we have used only the nuclear region. We use the lOGyr, 
lGyr, 100+10 Myr and 3Myr+FC as old, intermediate, young and FC 
components, respectively. 
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Figure 8. Same as Fig. [4] 



• ESO428-G014: CF04 carried out a study of the SP, in the 

spectral region between 3500 and 5200 A of the inner 200 pc of this 
source. Their study also used STARLIGHT. They found x y =20%, 
xi=47% and xo=i'i%. As can be observed in Tab[3]our NIR results 
agree with those obtained by CF04 in the optical. 

• Mrk 1066: The light in the optical region of this galaxy 
is dominated by young to intermediate age SPs, and the nu- 
clear spectrum is strongly dil uted by a non-thermal component 



clear spectrum is strong 
dGonzalez-Delgado et al.l 



ported by iRaimann et al 



20011). These results agree with those re- 



d2003l) . Our results show that light at 
1.22^im is dominated by a intermediate age SP diluted by a FC 



component, in agreement with the optical results (see Tab. [3]>- 
With respect to the stellar component, our results agree with 
iRamos Almeida et alj J2009T) . wich studied the same spectral re- 
gion by modeling the continuum with combinations of stars and 
black-body (with T=1000K) dilution, without the FC component. 

• Mrk 121 0: The optical SP of t his source has been studied by 
many authors dSchmitt et alj|l99^ : IStorchi-Bergmann et alj|200d, 
CF04) however, the results are controversial. While ISchmitt et all 
d 19991) argue that the opti cal light is dominated by a 10 G yr pop- 
ulation (54% at A 5870A). IStorchi-Bergmann et al.l J2000t> suggest 
that a young SP (or a power law) contributes with 50% of the flux 
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Figure 9. Same as Fig. [4] 



observed in A4020A. When modelling the spectral interval between 
A3500-5200 A, CF04 inferred that 53% of the flux observed in the 
spectrum Mrk 1210 was due to a young stellar population. They 
also found a contribution of 39% for the old component and 5% 
for the intermediate one. Our results disagree with those reported 
in the optical, as we found a lower fraction for the old component 
(~29%) and ~49% for the intermediate one. Moreover, ~20% of 
the flux belongs to young+fC components. 

• NGC 3227: We found a dominant young component (~ 40%) 
for th is galaxy. Our results agrees with those of lDavies et al.l d2007t 
l2006h which analyze the star formation in the inner 10 pc of this 



galaxy using the near infrared adaptive optics integral field spec- 
trograph SINFONI. They found that the light at this region is dom- 
inated by a 40Myr SP. 



• Mrk 573: The optical SP of this source was studied by 
ISchmittetaTI ( ll999h . They fit the W\ and continuum ratios, find- 
ing that 82% of the flux observed in A 5870A is due to a 10 
Gyr stellar popula tion. These results are in good agreement with 
those obtained by iRaimann et alj ( 120031) . IStorchi-Bergmann et al .1 
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Figure 10. Three continuum components of NGC 7674 and Mrk 124. Dot-short dashed line represents the SP (x y + x\ + x a ). The FC and hot dust component 
are represented by the long and short dashed lines, respectively. The solid line is the sum of the three components and the dotted line represents the observed 
spectrum. 
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Figure 11. Histograms comparing the population vector components. 
Light-weighted at the left panel and mass-weighted to the right one. 
The lines (colours) indicating each component of v are on the labels. 
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Figure 12. Histograms comparing our results in terms of light and 
mass fractions. At left we show a histogram of light- (filled) and mass- 
weighted (solid line) mean ages. At the right we present a histogram 
comparing the light (filled) and mass-weighted mean metallicities. 



( feOOOl) and lGonzalez -Delgado et al. I d200lh in the optical region^ 
Our NIR spectral synthesis points to a dominant intermediate age 
population (~53%) diluted by a FC component, which contributes 
with ~22 % of the observed cont i nuum flux. As for Mrk 1066, we 
agree with lRamos Almeida et al. with respect to the domi- 

nant presence of an intermediate age population, but we only found 
dilution by the FC component. 



• NGC 1097: Th e SP of this galaxy was studied in the UV by 
iBonatto et ail dl99oj) . They infer that the light at A2646 A is due to 
a fraction of ~40% of a young 16% to intermediate age and 

~44% to old SR These results agree with the detection of a young 
starbu rst in th e inner 9 pc of this galaxy l IStorchi-Bergmann et al.l 
I2005L see also lDavies etai] d2007l (20091) ) and with the optical SP 
studied by CF04, who found a contribution of 31% for the young 



9 These two latter groups have analysed the absorption lines species 10 We call young SP the sum of the contribut ions with age ^0.2Gyr. For 

present in the observed spectra. more details see Tab 1 1 of lBonatto et"aT] ll998l) . 
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component. However, there are some differences between the UV 
and optical SP for the intermediate age (CF found 37%) and for the 
old (CF found 12%, they also determine 19% for the FC) compo- 
nents. The results for the young SP of both groups of authors are 
consistent with our NIR synthesis. Regarding the intermediate and 
old age components, our results are more consistent with those of 
the optical region. The discrepancies between our results and those 
obtained in the UV are probably related to the fact that in the NIR 
we are integrating light through a deeper line of sight and thus, we 
detect old stars located more internally in the bulge of the galaxy. 

• NGC2110: iGonzalez-Delgado et alj feOOll) argue that the 
stellar optical absorption lines are similar to those of an old popula- 
tion. By fitting the nuclear continuum with an off-nuclear spectrum 
of the galaxy bulge as template, they conclude that the contribu- 
tion of a power law is unnecessary. In the case of an elliptical as a 
template spectrum, they needed to include a contribution ^25% of 
a power law (F\ cx A -1 ). These results are further confirmed by 
CF04, who found a contribution of 67% for the old component, 7% 
for the intermediate a ge and 26% for the young SP, as well as by 
iRaimann et alj J2003h who determine contributions of, 53%, 32% 
and 10% for the old, intermediate and young age components, re- 
spectively. Our results indicate a dominant intermediate age stellar 
population and a strong FC component. The detection of the latter 
component may be related to the fact that the NIR is less affected by 
dust extinction and is consistent with the detection of broad com- 
ponents in the emission lines, (see Sec. |5.2t . 

• NGC 5728: CF04 fitted the SP of this object and concluded 
that its light between 3500 and 5200 A is dominated by intermedi- 
ate age stars (51%), with contributions of 20% and 28% of young 
and old SPs, respectively. These fractions agree with our NIR anal- 
ysis, where we found a dominant contribution of the intermediate 
age SP (~70%) and ~30% for the old component. However, we did 
not detect the young component in our fitting process. This can be 
associate d to the fact that NIR light is dominated by intermediate 
age stars dMarastorj|2005h or to our spectral coverage that misses 
the calcium triplet abso rption, lines which are more sensitive to 
young SPs (see Fig. 4 of lRiffel et al]2008bl) . 

• NGC 5929: According to lGonzalez-Delgado et all d200ll) the 
nucleus of this source is dominated by old SP, a conclusion ob- 
tained by studying the W\ of optical absorption lines. They also 
conclude that if the off-nuclear spectrum is used as template to fit 
the nuclear spectrum of NGC 5929, n o signs of dilution ar e ob- 
served. Simmilar results are obtained bv lRaimann et al.1 d2003h who 
found contributions of, 61%, 27%, 8% and 4% for the old, interme- 
diate and young and FC components, respectively. These results 
fully match our NIR synthesis (FC=0%, x v =l.l%, z/=24.5% and 
xo=69.7%) 

• NGC 5953: One of the galaxies studied by CF04. They de- 
rived a dominant intermediate age optical SP (xi=1A%) with small 
fractions of young and old stars (x y =l% and xo=l%) and a con- 
tribution of 18% of the FC. The NIR SP of this galaxy is divided 
into two components, the old, dominant, (~65%) and the interme- 
diate age (~34%). The FC component does not contribute to the 
NIR flux. 

• NGC 7682: The optical light of this source is dominated by 
an old SP (86%) with no contribution of intermediate age stars 
(CF04). A small contribution of young populations (8%) and 6% 
of a FC component is reported by CF04. Our NIR results disagree 
with those of the optical: the SP is shared between old (~63%) and 
intermediate (~35%) age populations. 
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Figure 13. Comparison of the population vectors obtained in the NIR (this 
work) and in the optical (CF04), for 7 objects in common. The symbols 
indicate the population vectors. The full line is the identity line, the dotted 
lines represent ±10% deviation from the identity. Open and filled symbols 
are the flux and mass fraction, respectively. 

Our NIR synthesis and that in the optical (CF04) have been 
analysed by the same method (the STARLIGHT code), allowing 
for a proper comparison of the objects in common to both stud- 
ies. Fig. [13] summarises this comparison. In general, our results 
do not agree with those in the optical. Part of the differences can 
be accounted for by the fact that the NIR is more suitable for the 
identification of old SPs and the detection of the uni que absorp- 
tions re lated to 1 Gyr-old SP featured by TP-AGB stars jRiffel et all 
l2008bl) . This hypothesis is appropriate for the case of NGC 5953, 
where the differences between our results and those of CF04 oc- 
cur in the xj and xo components. In addition, the differences be- 
tween optical and NIR SPs can be associated to the fact that the 
NIR probes SPs buried deep in the dust. An example of this sit- 
uation is NGC 21 10, where we have detected the presence of hot 
dust in its integrated spectrum (Tab, [pi. For the case of NGC 1097, 
IStorchi-Bergmann et al.1 d2005l) report an obscuration of a central 
starburst, which they associate with a dusty absorbing medium. 

Our finding that the central regions of the galaxies contain a 
substantial fraction of intermediate-age SPs (see Figs. 1 1 1 landl 14t. 
together with the prolonged SFH, is very similar to the picture 
draw n in the case of central star-forming rings, based on optical 
data dAllard etal .1120061 : Is arzi et al. 1 120071) . This might support the 
scena rio where central s tar formation often occurs in circumnuclear 
rings dSarzi et al1l2005l : lshields et alj2007l) . 

5.1.2 The CN versus intermediate age stellar population 

Our NIR approach offers a unique opportunity to investigate in 
a more consistent way the relation between the CN molecular 
band and the unambiguous evide nce of an intermediate age SP 
( lMarastorj2005l ; |Riffel et al.l2007l ; ?). Fig.[T4lpresents a histogram 
comparing the i ntermediate age SP of the galaxies where CN was 
clearly detected dRiffel et"al]|2007|) and those with no detections in 
a visual inspection. The objects with a clear CN detection have con- 
tributions of the intermediate age component higher than ~20%, 
with a mean value of 40±15%. Regarding the galaxies where CN 
was not detected, the values are biased to contributions lower than 
~30%. The mean value derived for the intermediate age for the 
sources without CN detection is 27±20%. 

In four galaxies where CN was not detected we obtained frac- 
tions of the xi component higher than ~35%. This result may ap- 
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Figure 14. Histograms comparing the intermediate age component of the 
galaxies with CN detection (empty histogram) and un-detection (shaded). 



pear contradictory, as the CN band is the stellar feature that most 
suitably traces the intermediate age component. For NGC 5728, for 
instance, we associate the 69% of the intermediate age component 
with our shorter spectral coverage (see Sec. 15.1. ft and the low S/N 
ratio. For ESO428-G014, Mrkl210 and MCG 5-13- 17, the CN 
band is totally filled by the Pa7 emission line (Fig 2 o f iRiffel et al.l 
l2006h . 

We conclude this section arguing that our NIR spectral fitting 
confirms that the detection of the CN band in the spectrum of a 
galaxy can be taken as an unambiguous evidence of the presence of 
an intermediate age SR 



5.2 The featureless component 



As discussed in Sec l4.ll it is very difficult to distinguish a reddened 
young strarburst from a F v oc u~ power law. However, this ef- 
fect is even harder in the optical, where the main difference between 
a 5 Myr SSP and a FC seen through an absorption A v ~ 2 — 3 mag 
by dust is the presence of the Balmer absorption lines and Balmer 
jum p in the blue side in the former (CF 04). 

Icid Fernandas & Terlevicrj l fT995T) predicted that a broad com- 
ponent in H/3 becomes distinguishable whenever the scattered 
FC contributes with ^20% to the optical continuum light. With 
our synthesis we can investigate this issue in the NIR. Six out 
of 15 Sy 2 galaxies of our sample display a FC contribution 
higher than >20% (Tab. [3}- Interestingly, a broad component in 
the Hydrogen lines is detected in t he spectra of the six sources : 



Mrk 1066 (e.g. lVeilleux et al. II 19971) . Mrk573 jNagao et alj|2004l). 
NGC 1275 (e.g. iHo et al.lll997l), NGC 21 10 (eXllleunanen et all 
l2003f). NGC 262 (e.g. [Miller & Goodrich] 1 199Ch and NGC 7674 



(e.g. lMiller & Goodrichll 19901) . 

For two of our galaxies a broad component is reported in the 
literature, and we do not detect strong contribution of the FC com- 
ponent. The first one is Mrk 1210, wh ere a broad component of Ha 
and H/3 is detected in po larised light JTran et alj 1 19921 ; Tran Il995h 
and in the NIR Hi lines dMazzalav & Rodriguez-Ardilall2007l) for 
which we find a ~12% contribution of the FC component and 
xy ~ 10%. Our NIR synthesis for this galaxy is consistent with that 
obtained by CF04 in the optical, but we tend to find a lower con- 
tribution of the young component. The second object is NGC 5953, 



where iGoncalves et alj i ll 9991) report the possible detection of a 
very weak broad Ha line, but we did not find any contributions 
of FC and xy- CF04 report a contribution of 18% for the non- 
thermal comp onent. No broad c omponents are detected in our NIR 
spectrum (see lRiffel et alj2006h . 

The minimum F C contribution predicted by 

ICid Fernandes & TerlevicH dl995h seems to be reflected in 
the NIR, as we clearly detect FC fractions higher than 20% 
in the Sy 2 objects with a broad component in the Hi lines. 
Therefore, our results reinforce their predictions. In addition, we 
detect FC >20% in almost all Sy 1 sourcej"! which is consistent 
with the above arguments. However in four of our Sy 1 galaxies 
(Mrk291, Arp 102B, MCG-5-13-17 and NGC 1097) we detect 
fractions of FC <20%. For the first three objects we associate 
this ambiguity to the absence of features suitable to properly fit 
the absorption spectra and to a poor S/N ratio in the blue side of 
the spectrum (see below). As discussed by CF04 and CF05, a high 
S/N ratio is required for an adequate detection of the different 
contributions to the integrated spectra. Interestingly, the three 
objects display the presence of hot dust, wich is in full agreement 
with the nature of the Sy 1 objects predicted by the unified 
model for AGN (see Sec [53). NGC 1097, ori ginally classified 



as LINER on the basis of its optical s p ectrum l Keell 19831) 



iyclai 
ill 1983 
12003 



reclassified by IStorchi-Bergmann et alj ( 1 19931 , 1 19971 . 120031) as a 
Sy 1 after observing broad Balmer emission lines, a featureless 
blue continuum and double peak profiles. Our spectral synthesis 
for NGC 1097 agrees with that obtained by CF04. Moreover, our 
value o f xy ~20% is consistent with the starburst nature of this 
galaxy dStorchi-Bergmann et alj|2005l) . 

It is worth to mention at this point that the Wot values, mea- 
sured for almost all object of our sample (see Sec.[3j, even if the 
dilution is considered, are consitent with values measured in nor- 
mal spiral galaxies (<15A lBica & Alloin|[l987h . However, in two 
cases, NGC 262 and NGC 21 10, if we account for the dilution we 
reach very high values for WcaT (~ 17 A and ~25A, respectively). 
One explanation is that the WcaT of these sources is overestimated 
due to the effects on the continuum and absorption lines by telluric 
features (see Tab. [T). Other posibility lies in the fact that the FC 
component is more sensitive to shorter wavelengths (A < 10000A 
see Fig. |10t. which is our region with the lower number of con- 
straint^ and the lower S/N (the border of the spectrum; A ^ 9 000 
A). Therefore, the FC component can be overestimated in these 
sources. 



5.3 Dust emission 

A rapid look at column 5 of Table [3] shows an excess of emis- 
sion over the stellar population and the featureless continuum for 
half of the AGN spectra. This excess (see Fig. I lot, characterised 
by a Planck distributio n, suggests emission from hot dust grains. 
Since iBarvainisJ dl987t) . evidence of the presence of dust near the 
sublimation t emperature has been observed in the central region 
of AGNs (e.g.lMarco & AlloirJjl998L[2000l;|Rodriguez-Ardila etail 
120051 ; iRodriguez-Ardila & Mazzalavl 120061: iRiffel et alj|2009h ."ln 
order to quantify the dust contribution, we show in Tab.|4]the indi- 
vidual contribution of each blackbody over the temperature interval 



11 With a mean value of ~20% if all Sy 1 sources are considered and 
32% if we exclude Mrk 291, Arpl02B, MCG-5-13-17 and NGC 1097. 

12 Due to the large number of emission lines which are masked out. 
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800 - 1400 K. Note that as the T=700K and T= 1 000 K compo- 
nents were not detected, we dropped them from Tab. [4] 

Only in five Sy2 galaxies we have detected hot dust, while 
a positive detection is reported for all Sy Is. The obvious rea- 
son for this is that we are dealing with hot dust in the K- 
band (T ~1 000 K), and in the case of Sy2s the dust is cooler 
(T ~600K), thus more visible in the L or M -bands. This hypoth- 
esis is supported by the fact that in two Sy 1 galaxies, Mrk 1239 
and Mrk 766, we can see directly the presence of hot dust in the K- 
band feodriguez -Ardila et al .l2005l : lRodriguez -Ardila & Mazzalavl 
120061) . There are two remarkable cases, Mrk 573 and Mrk 1066, 
where we have detected a significant fraction of the FC compo- 
nent and no hot dust. There are some possible explanations for this, 
the first two are that discussed above. Another possibility is that at 
least a significant fraction of the detected FC component is due to 
a very young starbust (t <5 Myr). The latter, as discussed above 
(see Sec. l5,2t lies in the fact that the FC can be overestimated due 
the small number of constraints in the shorter wavelengths. 

The two main constituents o f interstellar du st are graphite 
(carbon) and silicate grains (e.g. iKriiegell |2003|) . As discussed 
in Sect. [5] the evaporation temperatures of gra phite and silicate 
grains are 1500 K and ~ 1000 K, respectively fearvainislll987l ; 
iGranato & Danesd 1 1994b . The temperature derived for the BB 
component of almost all objects is ^900 K (see BB C and BBh in 
Tab. [3), suggesting that hot dust close to the central source is prob- 
ably composed by graphite grains instead of silicates. The only dis- 
crepant objects are Mrk 1210 and NGC 1275 with T=800K. Con- 
siderin g, however, that th e spatial resolution in those objects is ~ 
400 pc jRiffel et"al]|2006b . and that the d ust temperature is a fu nc- 
tion of the distance to the central source dMarco & Allomlll998b . it 
is very likely that dust at higher temperatures exists closer to the 
central source. This would rule out the possibility of silicates as the 
main component of the nuclear dust grains. 

Assuming that the temperature of the different BB, found for 
each galaxy represents the distribution of temperatures for graphite 
grains in the nuclear region, we can estimate the lower limit of 
the hot dust mass responsible for the observe d if -band excess (see 
Fig. Hot following the approach developed bv lBarvaini sl dl987b . 

The inf rared spectral lu minosity of each dust grain can be ob- 
tained from jBar vainis 19§3): 

L a J >ir = 4n 2 a 2 Q„B„(T gr ) [ergs" 1 Hz" 1 ], (8) 

where a is the grain radius; Q v = qirV 1 is its absorption effi- 
ciency and B v (T gr ) is the Planck function for a grain at temper- 
ature Tg r . L g J ir was calculated assuming a typical grain radius 
a = 0.05u m dBM-yainisll 19871 ; iKishimotolkoOTb g lr =1.4x 10~ 24 
and 7=1.6 (Barvainis 1987|). The values of L 9 J ir are shown in 
Tab. |1 

The hot dust mass, Mh d, can be obtained by the equation 
dRodnguez-Ardila et alj2005b : 



Table 5. Hot dust masses found in AGNs. 



47T 3 

—a NhdPsi, 



(9) 



where Nhd 



is the number of hot dust grains and p gl - 

- HD 



is the density of the grain. L ir is the total NIR luminosity due 
to hot dust. It can be derived from the integrated the flux of each 
BB contribution over the spectral range between 0.01 and 160/im 
found on each galax\F 3 1 Then, we multiplied the integrated nor- 



Galaxy 


M HD (Mq) Spectr 


NGC 7582 


2.8 x 10" 


3 


yes 


Mrk 1239 


2.7 x 10" 


2 


yes 


Mrk 766 


2.1 x 10" 


3 


yes 


NGC 1068 


1.1 x 10" 


3 


no 


NGC 7469 


5.2 x 10" 


2 


no 


NGC 4593 


5.0 x 10" 


4 


no 


NGC 3783 


2.5 x 10" 


3 


no 


NGC 1566 


7.0 x 10- 


4 


no 


Fairall 9 


2.0 x 10" 


2 


no 



Riffel et al. (2009) 

Rodrfguez-Ardila & Mazzalay (2006) 

Rodriguez-Ardila et al. (20051 

Marco & Alloin (20001 

Marco & Alloin (19981 

Santos-Lleo et al. (19951 

Glass (19921 

Baribaud etal. (19921 

Clavel. Wamsteker & Glass (19891 



malized flux by the actual flux at 1.2/im (our normalisation point) 
and convert it to the adequate units (from erg cm" 2 s" 1 A" 1 to 
erg cm" 2 s" 1 Hz" 1 ). The final result of this process is the flux of 
each BB contribution (Fbb)- The values derived for our sample 
are presented in Tab. [4] The L^ r D was estimated using Fbb and cz 
listed in Tab.|4](we have adopted Ho = 75 kms - Mpc" 1 ). 

Finally, we have estimated the lower limi t for hot dust mass, 
for g raphite grains with p gI = 2.26 gem -3 dGranato & Danesel 



1 1994b and using Eq.|9] The hot dust mass of each BB distribution, 
as well as the total hot dust mass (^Mhd ), are presented in Tab.|4] 

To compare the mass values derived for our sample and those 
reported by other authors, Tab.[5]lists the ma sses available in litera- 
ture, determined following dBarvainisll 1987b . It is important to note 
that only three objects, Mrk 1239, Mrk 766 and NGC 7582, of the 
nine listed in Tab.[5]have the hot dust mass estimated by means of 
spectroscopy. In the remaining six, the masses were determined us- 
ing photometry. The spectroscopic approach allows a careful sub- 
traction of the power-law contribution and the stellar population. 
Moreover, our results have the advantage, over previous determina- 
tions, of considering the SP, the FC and the BB in the same fitting 
process. Thus, our study has increased significantly the number of 
AGNs with the mass of hot dust estimated. 

Given that the radius of the integrated region covered by most 
of our data is less than ~500 pj 1 " 4 ! this value sets an upper limit to 
the size of the hot dust emission region. We may further constrain 
this emission region if we consider that in half of our sample the 
spatial resolution is < 400 pc. It means that the bulk of the hot dust 
is more likely concentrated close to the central source. 

The origin of the hot du st can be the putative torus required by 
the un ified model for AGNs dAntonucci & Millerlll985l ; lAntonuccil 
1 1993b . which is a natural dust reservoir. This hypothesis is fur- 
ther supported by the detection of hot d ust in the i nner ~ 25 pc 
of NGC 7582 dRiffel et al.l [20091) and by Ijaffe et al] d2004b . who 
studied the Mid-infrared spectrum of NGC 1068, using interferom- 
etry with the Very Large Telescope Interferometer-WWI. The latter 
work shows that the 10 /im emission is due to hot dust at 800 K con- 
centrated in the central pc of this object, surrounded by cooler dust 
(T=300K) in scales of ~2-4pc. However, the detailed discussion 
of dust distribution around the central region of AGNs is beyond 
the scope of this paper and is left for a forthcoming work (Riffel et 
al. 2009, in preparation). 



J Note that we have weighted the distribution according to its contribu- 
tion to the total SED. 



Except for Arp 102 B, Mrk 124 and Mrk 291. 
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Table 4. Dust properties for our galaxy sample. 



Seyfert 2 Seyfert 1 



T(K) L gr . (a) 




NGC262 NGC2110 Mrkl210 Mrkl275 NGC7674 


Mrk334 MCG 5-13-17 Mrkl24 NGC3227 NGC4051 Mrk291 Arpl02B Mrk896 




Cont.(%) 








0.5 


0.3 































Fbb (b) 






7.53 


10.79 




















SOU 15.0s 


Lbb (c) 






2.62 


6.36 






















N 9 r (d) 






173.76 


421.51 






















M_s_D (e) 






103.33 


250.65 






















Cont.(%) 


1.8 











2.2 




























Fbb (b) 


10.51 








12.17 


















nnn i c m 


Lbs (c) 


4.54 








19.46 




















N ffr (d) 


175.21 








750.60 




















M/fo (e) 


104.19 








446.34 




















Cont.(%) 





























2.3 













Fbb (b) 




















7.12 








1 inn m i n 

1200 97.39 


Lbb (c) 




















0.07 










N gr (d) 




















0.76 










Mjfjy (e) 




















0.45 










Cont.(%) 














5.7 














2.0 













Fbb (b) 










2.54 










4.05 








1300 140.75 


Lbb (c) 










4.07 










0.04 










N 9r (d) 










28.91 










0.30 










Mac (e) 










17.19 










0.18 










Cont.(%) 


3.6 


0.9 











8.2 


3.9 


22.6 


2.1 


2.5 


4.5 


3.0 


13.7 




Fbb (b) 


1.21 


1.07 








5.17 


3.70 


6.54 


4.35 


3.61 


0.56 


1.13 


6.29 


1400 197.92 


Lbb (c) 


0.52 


0.12 








4.76 


1.10 


39.64 


0.12 


0.04 


1.32 


1.26 


8.39 




N 9r (d) 


2.64 


0.63 








24.04 


5.54 


200.31 


0.63 


0.19 


6.66 


6.36 


42.42 




M H d (e) 


1.57 


0.37 








14.29 


3.29 


119.11 


0.37 


0.11 


3.96 


3.78 


25.22 


EMflD (e) 




106 


0.37 


103 


251 


463 


14 


3 


119 


0.4 


0.75 


4 


4 


25 


cz(f) 




4507 


2335 


4046 


4046 


8671 


6579 


3731 


16878 


1157 


700 


10552 


700 


7922 



Table notes: (a) 1 X lO^crgs" 1 Hz -1 ; (b) 1 X 10" 27 ergcm- 2 s" 1 Hz" 1 ; (c) 1 X 10" 27 crgs- 1 Hz" 1 ; (d) 1 X 10 43 ; (e) 1 X 1O" 5 M0; (f) From 
NED. Note that as we have not detected contributions for the T=700 K, 2.64 and T=l 100 K BB components we left them out from the table. 



6 FINAL REMARKS 

In this work we investigate the NIR spectra of 24 Seyfert galax- 
ies (9 Sy 1 and 15 Sy2) observed with the IRTF SpeX, obtained 
in the short cross-dispersed mode. Our main focus was the stel- 
lar population, AGN featureless continuum and dust contribution 
properties, along the full wavelength coverage (0.8/im - 2.4/im). 
We have analysed the absorption features located in the NIR. The 
approach followed here is based on the STARLIGHT code, which 
considers the whole observed spectrum, continuum and absorption 
features. This is the first instance where STARLIGHT is applied to 
this wavelength range. Besides, in this work we also consider for 
the first time hot dust as an additional element base. 

The main results can be summarised as follows. We found evi- 
dence of correlation among the W\ of Si 1 1 .59/im x Mg 1 1 .58/im, 
equally for both kinds of activity. Part of the WNai2.2if*m and 
Wco 2.3 M m strengths and the correlation between W^Nai2.2i M m and 
Wms i i.58/im appears to be accounted for by galaxy inclination. 
For the 7 objects in common with previous optical studies (based 
on the same method of analyses), the NIR stellar population syn- 
thesis does not reproduce well the optical results. Our synthesis 
shows significant differences between Sy 1 and Sy 2 galaxies. The 
hot dust component is required to fit the A'-band spectra of ~80% 
of the Sy 1 galaxies, and only of ~40% of the Sy 2. Besides, about 
50% of the Sy 2 galaxies require a featureless component contribu- 
tion in excess of 20%, while this fraction increases to about 60% 



in the Sy 1. Also, in about 50% of the Sy 2, the combined FC and 
Xy components contribute with more than 20%, while this occurs 
in 90% of the Sy 1. This suggests recent star formation (CF05) in 
the central region of our galaxy sample. We found that the light at 
1.223/im in central regions of the galaxies studied here contain a 
substantial fraction of intermediate-age SPs with a mean metallic- 
ity near solar. Moreover, our analysis confirms that the 1 . 1 /im CN 
band can be taken as an unambiguous tracer of intermediate-age 
stellar populations. 

One consequence of this work - especially because of the si- 
multaneous fitting of SP, FC and hot dust components, allowing a 
proper analysis of each one of them - is a 150% (400% if only spec- 
troscopic studies are considered) increase in the number of AGNs 
with hot dust detected and the mass estimated. 

What emerges from this work is that the NIR may be taken as 
an excellent window to study the stellar population of Sy 1 galaxies, 
as opposed to the usually heavily attenuated optical range. Our ap- 
proach opens a new way to investigate and quantify the contribution 
of the three most important NIR continuum components observed 
in AGNs. 
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